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Abstract.  Time  series  analysis  was  used  to  develop  a  model  that  de- 
picts pulpwood  use  in  the  Northeast.  The  model  is  useful  in  forecast- 
ing future  pulpwood  requirements  (short  term)  or  monitoring 
pulpwood-use  activity  in  relation  to  past  use  patterns.  The  model 
predicted  a  downturn  in  use  during  1980. 


In  many  decisionmaking  situations,  past 
values  may  be  the  best  variables  to  use  when 
forecasting  future  values.  Developing  sto- 
chastic models  can  be  both  a  necessary  and 
difficult  endeavor.  Such  models  are  nec- 
essary in  an  industrial  society  because  we  live 
in  a  world  of  increasing  population  and 
limited  natural  resources  and  raw  materials. 
For  example,  to  satisfy  regional  require- 
ments for  paper  and  allied  products,  ade- 
quate supplies  of  pulpwood  and  woodpulp 
must  be  forthcoming,  or  substitute  materials 
must  be  found.  With  forecasts  of  expected 
pulpwood  demand,  industrial  decision- 
makers can  devise  strategies  for  obtaining 
the  pulptimber  resources  that  will  be 
needed. 

The  difficult  aspect  of  developing  a  pre- 
dictive model  is  that  variables  that  are 
currently  known  to  be  highly  related  to  pulp- 
wood use  are  assumed  to  be  important 


prognosticators  of  future  use.  In  addition  to 
the  forecaster's  inability  to  incorporate 
technological  and  other  change  into  his 
predictive  model,  he  is  also  vulnerable  to 
the  mathematical  constraints  of  his  model. 
Simple  time  series  models  may  only  be  able 
to  predict  broad  trends  within  the  current 
data  base,  while  other  types  of  models  wdth 
more  parameters  may  account  for  every 
aberration.  When  these  models  are  used  to 
project  as  many  as  ten  periods  into  the 
future,  the  simple  model  often  proves  to  be 
the  most  helpful  because  the  decisionmaker 
understands  its  limitations  and  is  able  to 
allow  for  possible  errors.  Although  the  time 
series  model  may  predict  actual  future 
values  closely,  it  can  only  forecast  future 
points  if  the  past  trends  and  cycles  are  re- 
peated. This  fact  points  to  another  important 
use  of  time  series  analysis;  the  technique 
allows  the  analyst  to  determine  whether 
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the  system  he  is  studying  is  operating 
normally.  By  using  his  model  as  a  standard 
for  comparison,  he  can  identify  abnormal 
periods. 

DATA  BASE 

For  many  years  the  American  Pulpwood 
Association  has  collected  monthly  statistics 
on  pulpwood  use  and  ending  inventory  from 
all  member  company  pulpmills  in  the  United 
States  and  has  published  monthly  sum- 
maries by  geographic  regions  (American 
Pulpwood  Assoc.  1977).  Pulpwood-use 
statistics  for  nine  states  of  the  Northeast 
region^  between  January  1967  and  August 
1979  were  used  for  developing  a  time  series 
model. 

A  plot  of  the  regional  pulpwood-use  sta- 
tistics revealed  that  there  is  a  predictable 
seasonal  pattern  (Fig.  1).  The  normal 
sequence  of  use  in  the  Northeast  has  been 
to  start  the  year  with  low-use  levels.  Pulp- 
wood use  proceeds  to  build  until  peaks  are 
reached  in  the  spring.  These  peaks  are  fol- 
lowed by  a  drop  in  use  as  pulpmill  workers 
go  on  their  summer  vacations  and  major  mill 

'Connecticut,  Maine,  Massachusetts,  New  Hampshire, 
New  Jersey,  New  York,  Pennsylvania,  Rhode  Island, 
and  Vermont. 


maintenance  is  done.  If  the  demand  for 
woodpulp  remains  high,  low  pulpwood  use 
in  July  is  followed  by  a  resurgence  in  use 
that  extends  into  the  fall.  As  winter  ap- 
proaches and  inclement  weather  hampers 
pulpmill  operations,  pulpwood  use  once 
more  declines. 

METHODS 

A  number  of  quantitative  forecasting  tech- 
niques are  available  for  making  predictions 
based  on  historical  analysis;  regression 
(Neter  and  Wasserman  1974),  the  U.S.  Cen- 
sus Seasonal  Adjustment  Program  (X-11) 
(U.S.  Bureau  of  Census  1967),  and  exponen- 
tial smoothing  (Brown  1959)  are  only  a  few. 
Recently  Box  and  Jenkins  (1970)  proposed 
a  structured  approach  to  time  series  model 
building  by  unifying  material  and  techniques 
that  have  been  available  for  a  long  time. 
In  general,  their  analysis  technique  accounts 
for  four  components  of  a  series,  namely, 
long-term  trends,  cyclical  effects,  seasonal 
effects,  and  random  variation.  The  first  three 
components  are  referred  to  as  the  "signal" 
and  the  random  variation  is  called  "noise." 
A  model  is  developed  which,  to  the  greatest 
extent  possible,  eliminates  the  signal  so 
that  only  the  noise  remains.  The  relative 


contribution  of  the  noise  is  determined  by 
the  autocorrelation  function.  Based  on 
standard  tests  such  as  Bartlett's  (1946)  ap- 
proximation or  the  Q  statistic  (Box  and 
Pierce  1970),  if  the  model  fits  the  observed 
data  except  for  random  variation,  then  it  is 
deemed  acceptable.  A  time  series  model, 
however,  can  only  be  validated  when  fore- 
casts are  compared  with  actual  future 
values. 

The  autoregressive  integrated  moving 
average  process  (ARIMA)  developed  by  Box 
and  Jenkins  is  a  powerful  tool  for  modeling 
time  series.  Contrary  to  other  model-build- 
ing techniques,  the  analyst  does  not  pick  a 
specific  model,  but  instead  eliminates  in- 
appropriate models  until  he  is  left  with  the 
most  suitable  one.  MINITAB  II  (Ryan  1978), 
a  general-purpose  statistical  computing 
system,  has  recently  added  a  computer 
algorithm  for  time  series  analysis  based  on 
ARIMA.  The  user  can  generate  a  series  of 
models  by  specifying  the  number  and  kind 
of  parameters.  The  program  provides 
sample  estimates  of  the  model  parameters 
with  their  standard  errors,  plots  of  auto- 
correlation functions  and  partial  autocor- 
relation functions,  and  forecasts  of  future 


observations  and  their  95  percent  con- 
fidence intervals. 

The  preferred  model  is  the  one  which  best 
represents  the  data  with  the  fewest  number 
of  parameters  that  are  consistant  with  the 
stochastic  structure.  Parameters  should  be 
statistically  significant  and  the  residuals 
should  not  be  serially  correlated. 

RESULTS 

The  model  that  best  depicts  monthly  pulp- 
wood  use  in  the  Northeast  is: 

(1-B12)  (l-0iB-02B2)Zt  =  (l-ei2Bi2)at 
where: 

B'"  =  backward  shift  operator  (B'"  Zt=  Zfm) 
Zt  =  observation  at  time  (t) 
at  =  random  error  at  time  (t) 

The  modeP  has  two  autoregressive  param- 
eters (01  and  02)  and  one  seasonal  moving- 


^For  detailed  explanation,  see  Box  and  Jenkins  1970, 
p.  8. 
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the  system  he  is  studying  is  operating 
normally.  By  using  his  model  as  a  standard 
for  comparison,  he  can  identify  abnormal 
periods. 

DATA  BASE 

For  many  years  the  American  Pulpwood 
Association  has  collected  monthly  statistics 
on  pulpwood  use  and  ending  inventory  from 
all  member  company  pulpmills  in  the  United 
States  and  has  published  monthly  sum- 
maries by  geographic  regions  (American 
Pulpwood  Assoc.  1977).  Pulpwood-use 
statistics  for  nine  states  of  the  Northeast 
region'  between  January  1967  and  August 
1979  were  used  for  developing  a  time  series 
model. 

A  plot  of  the  regional  pulpwood-use  sta- 
tistics revealed  that  there  is  a  predictable 
seasonal  pattern  (Fig.  1).  The  normal 
sequence  of  use  in  the  Northeast  has  been 
to  start  the  year  with  low-use  levels.  Pulp- 
wood use  proceeds  to  build  until  peaks  are 
reached  in  the  spring.  These  peaks  are  fol- 
lowed by  a  drop  in  use  as  pulpmill  workers 
go  on  their  summer  vacations  and  major  mill 
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maintenance  is  done.  If  the  demand  for 
woodpulp  remains  high,  low  pulpwood  use 
in  July  is  followed  by  a  resurgence  in  use 
that  extends  into  the  fall.  As  winter  ap- 
proaches and  inclement  weather  hampers 
pulpmill  operations,  pulpwood  use  once 
more  declines. 

METHODS 

A  number  of  quantitative  forecasting  tech- 
niques are  available  for  making  predictions 
based  on  historical  analysis;  regression 
(Neter  and  Wasserman  1974),  the  U.S.  Cen- 
sus Seasonal  Adjustment  Program  (X-11) 
(U.S.  Bureau  of  Census  1967),  and  exponen- 
tial smoothing  (Brown  1959)  are  only  a  few. 
Recently  Box  and  Jenkins  (1970)  proposed 
a  structured  approach  to  time  series  model 
building  by  unifying  material  and  techniques 
that  have  been  available  for  a  long  time. 
In  general,  their  analysis  technique  accounts 
for  four  components  of  a  series,  namely, 
long-term  trends,  cyclical  effects,  seasonal 
effects,  and  random  variation.  The  first  three 
components  are  referred  to  as  the  "signal" 
and  the  random  variation  is  called  "noise." 
A  model  is  developed  which,  to  the  greatest 
extent  possible,  eliminates  the  signal  so 
that  only  the  noise  remains.  The  relative 


contribution  of  the  noise  is  determined  by 
the  autocorrelation  function.  Based  on 
standard  tests  such  as  Bartlett's  (1946)  ap- 
proximation or  the  Q  statistic  (Box  and 
Pierce  1970),  if  the  model  fits  the  observed 
data  except  for  random  variation,  then  it  is 
deemed  acceptable.  A  time  series  model, 
however,  can  only  be  validated  when  fore- 
casts are  compared  with  actual  future 
values. 

The  autoregressive  integrated  moving 
average  process  (ARIMA)  developed  by  Box 
and  Jenkins  is  a  powerful  tool  for  modeling 
time  series.  Conti-ary  to  other  model-build- 
ing techniques,  the  analyst  does  not  pick  a 
specific  model,  but  instead  eliminates  in- 
appropriate models  until  he  is  left  with  the 
most  suitable  one.  MINITAB  II  (Ryan  1978), 
a  general-purpose  statistical  computing 
system,  has  recently  added  a  computer 
algorithm  for  time  series  analysis  based  on 
ARIMA.  The  user  can  generate  a  series  of 
models  by  specifying  the  number  and  kind 
of  parameters.  The  progi-am  provides 
sample  estimates  of  the  model  parameters 
with  their  standard  errors,  plots  of  auto- 
correlation functions  and  partial  autocor- 
relation functions,  and  forecasts  of  future 


observations  and  their  95  percent  con- 
fidence intervals. 

The  preferred  model  is  the  one  which  best 
represents  the  data  with  the  fewest  number 
of  parameters  that  are  consistant  with  the 
stochastic  structure.  Parameters  should  be 
statistically  significant  and  the  residuals 
should  not  be  serially  correlated. 

RESULTS 

The  model  that  best  depicts  monthly  pulp- 
wood use  in  the  Northeast  is: 

(1-B12)  (l-0,B-0,B-)Zt  =  (l-eivBi2)at 
where: 

B"'  =  backward  shift  operator  (B"'  Zj  =  Zt-m ) 

Zt  =  observation  at  time  (t) 

at  =  random  error  at  time  (t) 

The  model-  has  two  autoregressive  param- 
eters (0,  and  02)  and  one  seasonal  moving- 


=For  detailed  explanation,  see  Box  and  Jenliins  1970, 
p.  8. 
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average  term  (9i2).  Differencing  was  re- 
quired for  the  seasonal  portion  of  the  model. 
The  estimate  of  parameters  and  their 
standard  errors  are: 


Parameter 


01 

02 
012 


Estimated  Standard 
value  error 


0.6085 
0.2617 
0.8210 


0.0835 
0.0853 
0.0750 


The  pulpwood-use  model  was  judged  satis- 
factory in  that  all  of  the  parameter  estimates 
were  statistically  significant,  and  none  of 
the  sample  autocorrelations  were  statis- 
tically significant  so  that  the  residuals  were 
not  autocorrelated. 

Based  on  this  model,  Figure  1  shows  a 
comparison  of  the  actual  pulpwood-use 
statistics  and  the  predicted  values  for  the 
Northeast  between  January  1968  and  August 
1979.  In  addition,  pulpwood  use  is  fore- 
cast to  August  1980.  During  the  11-year 
period,  all  of  the  predicted  values  were 
within  ±13  percent  of  the  actual  values,  ex- 
cept in  July  1978  when  actual  pulpwood 
use  dropped  25  percent  below  predicted 
use.  This  illustrates  how  the  model  can  be 
used  as  a  means  of  identifying  abnormal 
levels  of  use. 

DISCUSSION 

By  using  the  ARIMA  technique,  a  time 
series  can  be  fitted  with  a  mathematical 
model  that  is  optimal  in  the  sense  that  it  as- 
signs smaller  errors  to  history  than  any  other 
model.  The  type  of  model  is  identified  and 
the  parameters  are  then  estimated.  This 
statistical  routine  is  currently  considered 
to  be  one  of  the  most  accurate  portrayals 
of  time  series  for  short-term  forecasting. 
Although  forecasts  are  used  by  resource  and 
raw '  material  planners  to  ebtimate  future 
req'^uirements,  they  are  also  used  by  econ- 
omists to  gauge  market  supplies  and 
demands  relative  to  the  norm.  Calendar  year 
1975  was  an  atypical  year  for  pulpwood 
use  in  the  Northeast.  Forecasts  based  on 
time  series  analysis  data  between  1967  and 


1974  tracked  monthly  use  of  between  450 
and  550  thousand  cords.  Actual  monthly 
use,  however,  dropped  to  between  400  and 
450  thousand  cords  during  1975.  As  a  re- 
sult, there  were  pulpwood  inventory  build- 
ups in  many  northeastern  woodyards.  Addi- 
tional wood  inventories  tie  up  working 
capital  and  increase  the  risk  of  wood  de- 
terioration during  the  longer  holding  period. 

Our  time  series  model  forecasts  that  7.0 
million  cords  of  pulpwood  will  be  required 
in  the  Northeast  between  August  1979  and 
August  1980.  This  figure  represents  a  1.0 
percent  decrease  in  use  from  the  same  period 
in  1978-1979. 
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